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Abstract 

Compact stars can have either hadronic matter or can have exotic states of matter hke strange 
quark matter or color superconducting matter. Stars also can have a quark core surrounded by 
hadronic matter, known as hybrid stars (HS). The HS is likely to have a mixed phase in between 
the hadron and quark phase. Observational results suggest huge surface magnetic field in certain 
neutron stars (NS) called magnetars. Here we study the effect of strong magnetic field on the 
respective EOS of matter under extreme conditions. We further study the hadron-quark phase 
transition in the interiors of NS giving rise to hybrid stars (HS) in presence of strong magnetic 
field. The hadronic matter EOS is described based on relativistic mean field theory and we include 
the effect of strong magnetic fields leading to Landau quantization of the charged particles. For 
the quark phase we use the simple MIT bag model. We assume density dependent bag pressure 
and magnetic field. The magnetic field strength increases going from the surface to the center of 
the star. We construct the intermediate mixed phase using Glendenning conjecture. The magnetic 
field softens the EOS of both the matter phases. The effect of magnetic field is insignificant unless 
the field strength is above lO^^G. A varying magnetic field, with surface field strength of lO^^G 
and the central field strength of the order of 10^''G has significant effect on both the stiffness and 
the mixed phase regime of the EOS. We finally study the mass-radius relationship for such type 
of mixed HS, calculating their maximum mass, and compare them with the recent observation of 
pulsar PSR J1614-2230, which is about 2 solar mass. The observations puts a severe constraint on 
the EOS of matter at extreme conditions. The maximum mass with our EOS can reach the limit 
set by the observation. 
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INTRODUCTION 



The central density of neutron stars exceed the nuclear saturation density (no ~ 
O.lSfm^^), thereby raising the idea that compact stars might contain deconfined and chirally 
restored quark matter in them. Recently, [l| the mass measurement of millisecond pulsar 
PSR J1614-2230 has set a new robust mass limit for compact stars to be M = 1.97±O.O4M0. 
This value, together with the mass of pulsar J1903+0327 of M = 1.667 ± 0.021 Mq Q is 
much larger than any of the highest precisely measured pulsar mass. These measurement 
has set for the first time a very strong limit on parameters of the EOS, which describes 



matter under extreme conditions 



After the discovery of pulsar [5] and connecting them with NS [6|, various EOS for nuclear 



a, 



matter has been proposed and refined |7H10|. The quark sector is not much well understood 



as the nature of strong interaction at extreme condition still remains a challenge. The 



strange quark matter (SQM) conjecture by Itoh, Witten UJ, |l2| consisting of almost equal 



number of up (u), down (d) and strange (s) quarks was supported by model calculations 



13| . The most simple and popular model which describes the properties of quark matter at 
such high densities is the MIT bag model [l^. New refined models based on results from 



recent experiments in laboratories has been proposed |15Nl7l|. Thus normal nuclear matter 
at high density and/or temperature is likely to be unstable against stable SQM and would 
eventually decay. 

Compact objects therefore can be made of either nuclear matter or quark matter. Stars 
which has only nuclear matter are called neutron stars (NS). Broadly there can be two classes 
of compact stars with quark matter. The first is the so-called (strange) quark stars (SS) 
of absolutely stable strange quark matter. The second are the so-called hybrid stars (HS), 
along with the hadronic matter they have quarks matter in their interior either in form of a 
pure strange quark matter core or color superconducting matter. In between the quark and 
the hadronic phase a quark-hadron mixed phase exists. The size of the core depends on the 
critical density for the quark-hadron phase transition and the EOS describing the matter 
phases. 

Usually, the presence of strangeness in quark and hadronic matter provides an additional 
degree of freedom and softens the EOS and therefore quark and hybrid stars cannot reach 
high masses. Thus the mass measurement of pulsar PSR J1614-2230 puts forward a strong 



constraint on such EOS. However, studies found that effects from the strong interaction, 
such as one-gluon exchange or color-sup er conduct ivit y can stiffen the quark matter EOS 



18 



23[. The first studies on the imphcations of the new 



2J] and 25|. They, however 



and increase their maximum mass 

mass hmits from PSR J1614-2230 for quark matter was done by 
did not include the effects from color-superconductivity. 

The presence of magnetic field in compact stars has an important role in astrophysics. 
New observations suggests that in some pulsars the surface magnetic field can be as high as 
10^^ — lO^^G. It has also been attributed that the observed giant flares, SGR 0526-66, SGR 



1900+14 and SGR 1806-20 



26|, are the manifestation of such strong surface magnetic in 



those stars. Such stars are separately assigned as magnetars. If we assume flux conservation 
from a progenitor star, we can expect the central magnetic field of such stars as high as 
10^'' — 10^®G. Such strong fields are bound to effect the NS properties. It can modify the 

or it can modify the EOS of matter of the star. The e : 
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- |34| and quark matter 
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metric describing the star 

of strong magnetic field, both for nuclear matter 
has been studied earlier in detail. 

Motivated by recent observations of maximum mass limits of compact stars and strong 
magnetic field in magnetars, in this work we want to explore their implications on the EOS 
of both phases of matter that may be present inside a neutron star. We study the hadron- 
quark phase transition inside a compact star with a mixed phase region in between the quark 
core and nuclear outer region. The paper is organized as follows. In Sec. II we discuss the 
nuclear EOS and the effect of Landau quantization due to magnetic field on the charged 
particles. In Sec. Ill we employ the simple MIT bag model for the quark matter EOS and 
the effect of magnetic field on the quarks (also due to Landau quantization). In Section IV 
we develop the mixed phase region by Glendenning construction. We show our results in 
section V for the density dependent bag constant and varying magnetic field for the mixed 
HS. Finally we summarize our results and draw some conclusion in section VI. 



MAGNETIC FIELD IN HADRONIC PHASE 

At normal nuclear density the degrees of freedom for the EOS are hadrons. To describe 
the hadronic phase, we use a non-linear version of the relativistic mean field (RMF) model 
with hyperons (TMl parametrization) which is widely used to construct EOS for NS. In this 
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model the baryons interact with mean meson fields [38|-|43|. 

The Lagrangian density including nucleons, baryon octet (A, S°'='=, and leptons is 
given by 

b ^ 

- (rrib - gabcr)]ipb + - {d^ad'^a - m^cr^) 

11 11 
+ Y.^L[il,d^-mL]iJL. (1) 

L 

Leptons C are treated as non-interacting and baryons h are coupled to the scalar meson a, 
the isoscalar- vector meson and the isovector- vector meson p^. There are five constants in 
the model that are fitted to the bulk properties of nuclear matter. This model is good enough 
to describe nuclear matter and the nuclear saturation point. But it is insufficient for the 
hyperonic matter, because the model does not reproduce the observed strong AA attraction. 
This defect can be remedied by adding two new meson fields with hidden strangeness, 
namely, the iso-scalar scalar a* and the iso-vector vector 0, which couple to hyperons only 



42|. 



The effective baryon mass is given by 

rrib* = mi,- g^a - g^*a*. (2) 
For the beta equilibrated matter the conditions is 

Pi = hpB + qiPe, (3) 

where hi and are the baryon number and charge (in terms of electron charge) of species z, 
respectively, /i^ is the baryon chemical potential and p^ is the electron chemical potential. 
For charge neutrality, the condition is 

Pc = Y.^i^i^ (4) 
i 

Hi is the baryon number density of particle i. 

The magnetic ffeld is assumed to be in the z direction, ^ = Bk. Now the motion of the 
charged particles are quantized in the perpendicular direction of the magnetic field. The 



landau quantized energy is given by [44 1 



Ei = Jpi^ + nii"^ + \qi\B{2n + s + l). 



(5) 



In the above equation n is the principle quantum number, s is the spin of the particle (either 
(+) or (-)) and pi is the momentum component along the field direction of particle i. We 
can write 2n + s + 1 = 2z/, where z/ = 0, 1, 2..., so that now the energy can be written as 



Pi^ + nii'^ + 2v\qi\B = Jpi^ + mf^ 



(6) 



where the i/ = state is singly degenerate. It should be remembered that for baryons the 
mass is rrib* . 

At zero temperature and in the presence of a constant magnetic field S, the number and 



energy densities of charged particles are given by 
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27r2 



and 



47r2 



E 



(7) 



(8) 



J, is the Fermi momentum for the level with the principal quantum number n and spin s 
and is given by 

p)l = E)'-ml\ (9) 

The Fermi energies are fixed by their respective chemical potentials. 

The number, energy, and scalar number densities of the neutral particles are given by 
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The total energy density of the system can be written as 
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(13) 



where the last term is the contribution from the magnetic field. The general expression for 
the pressure is given by 

P = J2iJ,ini- e. (14) 

i 

At the outermost surface of the star, that is at lower densities, the matter is composed of 
only neutrons, protons and electrons. Hence, at the low density regime, only the electrons 
and protons are affected by the magnetic field. Electron being highly relativistic, the number 
of occupied Landau levels by electrons is very large. The field strength under consideration 
is larger than the critical field strength of electron by several orders but very less than 
the critical field strength of protons. Therefore, the number of occupied Landau levels by 
protons is large. As the magnetic field increases with the increase of density, the number of 
occupied Landau levels gradually decreases for every species. 



MAGNETIC FIELD IN QUARK PHASE 

Considering the simple MIT bag model for the quark matter in presence of magnetic field 
we assume that the quarks are non-interacting. The current masses of u and d quarks are 
extremely small, e.g., 5 and 10 MeV respectively, whereas, for s-quark the current quark 
mass is taken to be 150 MeV, unless otherwise stated. 

For the same constant magnetic field configuration along the z-axis, the single energy 



eigenvalue is given by 44 1 



Ei = ^p^ + m,^ + 2u\q,\B (15) 
Then the thermodynamic potential in presence of strong magnetic field B{> B^^\ critical 



value discussed later) is given by [45 1 



For the zero temperature, the Fermi distribution is approximated by a step function and by 
interchanging the order of the summation over u and integration over E, 



a = E / , , , , dE^^El -m\- 2u\q,\B. (17) 

The upper limit of v sum can be obtained from the following relation 

Pf/ = -rrii^ -2v\qi\B>Q, (18) 



where pf^i is the Fermi momentum of the particle i. 

The upper hmit is not necessarily same for all the components. For a certain critical 
magnetic field strength the energy of a charged particle changes significantly in the quantum 
limit. For an electron with mass 0.5 MeV, the critical field strength is ~ 4.4 x 
whereas for a light quark (u or d), this value becomes ~ 4.4 x lO^^G, and for s-quark of 
current mass 150 MeV, it is ~ lO^^G. A compact star becomes unstable if the magnetic field 
strength becomes much greater than ~ lO^^G, and so many authors have neglected quantum 



mechanical effect of magnetic field on s-quarks 46| but in our calculation we include the 
quantum mechanical effect for all particles. 

Assuming the strange quark matter also to be charge neutral and in chemical equilibrium, 
we may write as 

/id = Ais = = /i« + A^e, (19) 
2n„ -Ud-ris- Sue = 0. (20) 
The baryon number density is given by 

nb = ^{nu + nd + ns). (21) 

Solving the above eqs f|T9l [20l [2T|) numerically, we obtain the chemical potentials of all the 
flavors and electron. Zero temperature approximation gives the number density of the species 

i {u, d, s, e) 

nj = TT 

47r2 



J2Jfi,^-m,^~2u\q,\B. (22) 



The total energy density and pressure of the strange quark matter is given by 

i i 

p = -J2n,-BG, (23) 

i 

where Bq is the bag constant. 

PHASE TRANSITION AND MIXED PHASE 



With t 
struction 



le above given hadronic and quark EOS, we now perform the Glendenning con- 



47( 1 for the mixed phase, which determines the range of baryon density where both 



phases coexist. Allowing both the hadron and quark phases to be separately charged, and 



still preserving the total charge neutrality as a whole in the mixed phase. Thus the matter 
can be treated as a two-component system, and can be parametrized by two chemical po- 
tentials, usually the pair (/ie,Ain), i-e., electron and baryon chemical potential. To maintain 
mechanical equilibrium, the pressure of the two phases arc equal. Satisfying the chemical 
and beta equilibrium the chemical potential of different species are connected to each other. 
The Gibbs condition for mechanical and chemical equilibrium at zero temperature between 
both phases is given by 

-PHp(/^e, /^n) = -PQp(/^e, l^n) = PmP- (24) 

This equation gives the equilibrium chemical potentials of the mixed phase corresponding 
to the intersection of the two phases. At lower densities below the mixed phase, the system 
is in the charge neutral hadronic phase, and for higher densities above the mixed phase the 
system is in the charge neutral quark phase. As the two surfaces intersect, one can calculate 
the charge densities pf^ and p^^ separately in the mixed phase. If % is the volume fraction 
occupied by quark matter in the mixed phase, we have 

XP^ + (1 - x)pT = 0- (25) 

Therefore the energy density cmp and the baryon density ^mf of the mixed phase can be 

obtained as 

cmp = XCqp (1 - x)eHP, (26) 
^MP = X^QP + {I- x)nnp- (27) 

RESULTS 

In neutron stars, the central part of the star has maximum density, therefore, it is much 
likely that the matter there undergoes a phase transition. As the density decreases towards 
the surface there is a probability of having nuclear matter and so in the intermediate stage 
there is a mixed phase, and as we go outwards we only have nuclear matter. The crust 
consisting mainly free electrons and nuclei which completes the star structure. 

The parametrization of the EOS of the hadron and quark phase is responsible for charac- 
terization of the mixed phase region. For the hadronic EOS we assume a fixed parameter set 
TMl, which reproduces the nuclear matter properties at high density quite well. However 
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FIG. 1. Pressure vs energy density plot with bag pressure of 170 and ISOMeV. 
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FIG. 2. Pressure vs number density plot with bag pressure of 170 and ISOMeV. 
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FIG. 3. Energy density vs number density plot with bag pressure of 170 and ISOMeV. 



the quark EOS can be controlled by changing the quark masses and the bag constant. The 
masses of the light quarks (u and d) are bounded and we take them to be 5 and lOMeV, re- 
spectively. The mass of s-quark is still not established, and can vary between 100 — SOOMeV, 
and we take them to be ISOMeV, unless otherwise stated. We regulate the bag constant 
(Bg) to characterize the mixed phase region. Now we use the Glendenning approach to 
construct the mixed phase, and obtain pressure vs energy density relation as given in fig [TJ 
In fig[T|, we have plotted the mixed phase EOS (pressure vs energy) with bag pressure 170 
and ISOMeV. For simplicity, we will denote Bq^^^ = 170MeV = Bg. The lower portion of 
the curve is nuclear phase (dotted line), the intermediate portion is the mixed phase (bold 
line) and the upper region is the quark phase (broken line). The curve with bag constant 
170MeV is much stiffer than the curve with bag pressure ISOMeV, because the bag pressure 
is negative to the matter pressure, making the effective pressure less. In fig|2]we have plot- 
ted pressure vs number density, and we find that the qualitative variation in the curves is 
same as that of fig[ll For bag constant 170MeV the mixed phase region starts at density 
0.2/m~^ and ends at 0.76/m~^. With bag constant ISOMeV the mixed phase region starts 
at density 0.22/m~^ and ends at O.S9/m~'^. In fig [3] we have plotted for the energy density 
vs number density, and we find a smooth curve, which does not differ from each other much. 
It is clear from the above figures that the main variation is due to the pressure, therefore, 
we only plot the pressure vs energy density curve. The above curves shows that as the bag 
pressure increases the range of mixed phase region increases, and there is a slight kink in 
the curve from going to the quark phase from the mixed phase. The EOS (or the pressure) 
for the nuclear matter is usually much stiffer than quark matter. As the bag constant with 
170 MeV is more stiffer than ISOMeV the kink in the former is much sharper than the latter 
one. Also as the latter curve is much flatter and so the mixed phase region is much extended 
there. By the Glendenning construction, we find that for a given mixed phase to exist the 
bag constant must be in between 170MeV and ISO MeV. 

The introduction of the magnetic field changes the EOS of the matter. The single par- 
ticle energy is now Landau quantized, and thereby it changes all the other thermodynamic 
variable of the EOS, namely the number density, pressure and the energy density. In fig HI 
we have plotted EOS for bag constant 170MeV with and without the magnetic field. The 
effect of magnetic field is insignificant when the field strength is less than 10^'^G, and also for 
this case the effect in the nuclear phase is very small. The magnetic field effect is less for the 

10 



500 
400 
> 300 

s 

3 200 
o. 

100 



"0 250 500 750 1000 1250 1500 

energy density (MeV/fm') 

FIG. 4. Pressure against energy density plot with bag pressure of 170MeV with and without 
magnetic field. The magnetic field is -B = 10^^ G. 
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FIG. 5. Pressure against energy density plot with bag pressure of ISOMeV with and without 
magnetic field. The magnetic field is i? = 10^^ G. 

nuclear matter than quark matter because, the nuclear EOS is much steeper than the quark 
EOS, thereby requiring much greater field to have any sound effect. We have plotted the 
same for bag constant ISOMeV (fig [5]), and for comparing the two bag constants, we have 
plotted figO For the bag constant 170 MeV, the EOS curve with magnetic field extends 
up to density 0.8/m~^, and for 180 MeV it extends upto density 0.92/m~^. The change 
in the mixed phase region is about 5 — 7%. Magnetic field makes the curve softer due to 
the negative effect of landau quantization on the matter pressure and the positive effect on 
the matter energy density. As shown in the figures, the effect of magnetic field begins to 
appear on the EOS of the matter when the field strength is above lO^^G. Such field has very 
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FIG. 6. Pressure against energy density plot with bag pressure of 170 and ISOMeV. For the bag 
pressure 170MeV we have also plotted curve with magnetic field of strength lO^^G. 

less effect on the nuclear matter but has considerable effect on the mixed and quark matter. 
With the onset of the magnetic field the mixed phase region gets extended. The magnetic 
field makes the quark matter EOS more flatter and therefore the mixed phase region is much 
extended. For fixed bag constant and fixed magnetic field value throughout, we cannot go 
to field strength above lO^^G, as it is bounded by observation of surface magnetic fields in 
magnet ars. 

Next we assume a density dependent bag constant. In the literature there are several 
attempts to understand the density dependence of bag constant 



48 



49|; but still there 



is no definite picture, and most of them are model dependent. We parametrized the bag 
constant in such a way that it attains a value -Boo, asymptotically at very high densities. The 
experimental range of Boo is given in Burgio et al. 50, |5lJ, and from there we choose the 
value = ISOMeV. With such assumptions we then construct a Gaussian parametrization 



given as 



50 



5l| 



Bgnirib) = Boo + {Bg - Boo) exp 



(28) 



The lowest value B^o, is the lowest value of bag pressure which it attains at asymptotic high 
density in quark matter, and is fixed at ISOMeV. The quoted value of bag pressure, is the 
value of the bag constant at the nuclear and mixed phase intersection point denoted by Bg in 
the equation. The value of Bgn decrease with increase in density and attain B^o = ISOMeV 
asymptotically, the rate of decrease of the bag pressure is governed by parameter 13. 

The observed magnetic field of the magnetars is of the order of ~ 10^^ — lO^^G. The flux 
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FIG. 7. Pressure as energy density plot with bag pressure of 170MeV, with and without varying 
magnetic field. The varying magnetic field has Bq = lO^^G and a = 0.005. 

conservation of the progenitor star may give the central field as high as ~ 10^^ — lO^^G.We 
assume that the parametrization of the magnetic field depends on the baryon number density. 



Therefore we assume a simple density dependence, given by 



30 



3l| 



(29) 



where a and 7 determines the magnetic field varaiation for fixed surface field Bg and asymp- 
totic central field Bo. The value of B depends mainly on Bq, and is quite independent of Bg. 
Therefore we vary Bq, whereas surface field strength is taken to be fixed at Bg = lO^^G. We 
keep 7 fixed at 2, and vary a for to have the field variation. Previous authors considered 
very high magnetic field value at the center, few times 10^*^0, but we would assume the 
maximum field to be of the order of few lO^^G. As this is somewhat low value from other 
previous assumptions, but it is more likely to be present in most magnetars. 

Next, we vary the magnetic field, and the bag constant is kept constant. First we consider 
the bag constant to be Bg = 170MeV. In fig [7] we have plotted curves for zero magnetic 
field and with Bq = lO^'^G with a = 0.005. As we vary the magnetic field, the magnetic 
field increases as we go towards to the center of the star. The field quoted in the figure is 
asymptotic field value. With Bq = and a = 0.005, the field strength is 4 x lO^^G at 

lOriQ. It is clear from the figure as the field strength increases, the curve becomes less stiffer. 
The change in the curve stiffness is due to the fact that the magnetic pressure due to landau 
quantization act in the opposite direction of the matter pressure, whereas, for the magnetic 
stress it acts towards the matter energy density. The two effect reduces the stiffness of 
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FIG. 8. Pressure against energy density plot with constant and varying bag pressure, Bg = 
ITOMeV. 
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FIG. 9. Pressure against energy density plot having density dependent bag pressure Bg = 170MeV, 
with and without magnetic field. The magnetic field strength is i? = 2 x lO^^G. 

the EOS (pressure vs energy density curve). It is also clear that the nuclear region (the 
low density regime) is not much affected by the magnetic field as there the magnetic field 
strength is low, whereas the quark sector (higher density regime) is the most effected region 
as the field strength is maximum there. However, the mixed phase region is moderately 
affected (the intermediate region). 

In fig [H] we plot curves with and without varying bag pressure, Bg = ITOMeV. For the 
curve with variation, at higher densities the bag pressure decreases, making the effective 
matter pressure higher. Therefore the pressure against energy density plot for this case is 
much stiffer. Also the mixed phase region gets shrunken due to the varying bag pressure. 
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FIG. 10. Pressure against energy density plot having density dependent bag pressure Bg = 170MeV 
and ISOMeV. Also shown in the figure the magnetic field {B = 2 x lO^^G) induced EOS curve for 
Bg = ITOMeV. 
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FIG. 11. Pressure with energy density plot having density dependent bag pressure ITOMeV, without 
magnetic field and with two different magnetic fields, having a = 0.005. 

The mixed phase region now only extends up to density 0.53/m~^. The change in the mixed 
phase region is about 40%. Therefore the change in the mixed phase region is much more 
influenced by varying bag pressure than due to magnetic field. We have plotted fig |9] with 
varying bag pressure Bg = ITOMeV, with and without constant magnetic field. The magnetic 
field employed for this plot is 2 x lO^^G. The change in the slope of the curves is due to the 
Landau quantization effect. The magnetic pressure acts opposite to the matter pressure, 
making the curve flat. For comparison, we have plotted flglTOlwith density dependent bag 
pressure, Bg = 170 and Bg = ISOMeV, and obtain quantitative same result. 
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FIG. 12. Pressure with energy density plot having density dependent bag pressure 170MeV, without 
magnetic field and with two different magnetic fields, having a = 0.01. 
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FIG. 13. Pressure with energy density plot having density dependent bag pressure 170MeV, without 
magnetic field and with same magnetic field but different a values. 

The curves for which both the bag constant and the magnetic field varies are of utmost 
importance. Fig [TT] shows curves for varying bag pressure 170MeV, without magnetic field 
and with varying magnetic field, Bq = lO^^G and 2 x lO^^G having a = 0.005. For the 
above values the field strength is 4 x lO^^G and 7.8 x lO^^G, at density IOuq. As the value 
of Bq increases the slope of the EOS curves becomes more and more soft, because the value 
of magnetic pressure increases with increase in field strength. As the magnetic pressure 
increases the effective pressure decreases making the curves flatter. In fig [12] we plot the 
same set of curve only for a = 0.01. With such a value, the asymptotic Bq = lO^^G gives 
field strength of 6 x lO^^G at lOno baryon density. For Bq = 2 x lO^^G the field strength 
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FIG. 14. Pressure vs energy density plot with two different density dependent bag pressure IGOMeV, 
ITOMeV. We fiave also plotted the magnetic field induced (field strength Bq = 2 x lO^^G) EOS 
curve for bag pressure IGOMeV having a = 0.01. 

is 1.21 X lO^^G at the same lOno baryon density. As the variation (a) becomes stiffer, the 
EOS curve becomes softer. This is seen clearly in fig [121 

We find for such varying bag constant and varying magnetic field, the change in the 
curves from the non varying non magnetic case is maximum. There is considerable change 
in the stiffness of the curves and also change in the mixed phase region. Towards the center, 
the magnetic field increases whereas the bag pressure decreases. On one hand the low bag 
pressure makes the curve stiffer whereas on the other hand large magnetic field strength 
makes the curve flatter. The low bag constant makes the mixed phase region to shrink, 
and the larger magnetic field tries to expand the mixed phase region. The effect of bag 
pressure is greater than the magnetic field and therefore the mixed phase is smaller than 
the constant bag pressure case. On the low density side, the effect of magnetic field is 
insignificant. Therefore the phase boundary between the nuclear and mixed phase is not 
much affected. 

For a varying bag constant we can have a significant mixed phase region with Bg = 
160MeV (fig fT4l) . The curve with bag pressure 160MeV is stiffer than other curves. This is 
because the bag pressure of Bg = 160MeV is lower than other higher bag pressure. Therefore, 
the effective matter pressure for this curve is higher than any other curve, which is reflected 
in the stiffness of the curve. For bag constant 160MeV the mixed phase region starts at 
density 0.15/m~^ and ends at 0.38/m~^. 
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FIG. 15. Gravitational mass (in solar mass) against radius plot of a star sequence with two different 
density dependent bag pressure, 170 and ISOMeV. 

Assuming the star is non rotating and has spherically symmetric, the distribution of 
mass is in hydrostatic equilibrium. The equilibrium configurations solution are obtained by 
solving the Tolman-Oppenheimer-Volkoff (TOV) equations js^l] for the pressure -P(e) and 
the enclosed mass m, 



dP{r) Gm{r)e{r) [1 + P(r)/e(r)] [1 + 47rr^P(r)/m(r)] 



dr 
dm{r) 
dr 



1 — 2Gm{r)/r 



Anr e(r) 



(30) 
(31) 



G is the gravitational constant. Starting with a central energy density e(r = 0) = ec, we 
integrate out until the pressure on the surface equals the one corresponding to the density 
of iron. This gives the stellar radius R and the total gravitational mass is then 

rR 



Mg = m{R) = Att dr r e{r) . 
Jo 



(32) 



For the description of the NS crust, we have added the hadronic equations of state with the 
ones by Negele and Vautherin 53| in the medium-density regime, and the ones by Feynman- 



Metropolis- Teller 



54l | and Baym-Pethick-Sutherland |55j for the outer crust 



Fig [TB] shows the gravitational mass M (in units of solar mass Mq) as a function of 
radius R, for varying bag pressure Bg = 170 and ISOMeV. As the bag pressure increases the 
curve becomes flat as the effective matter pressure decreases (bag pressure being negative) 
thereby decreasing the maximum mass of the star. We notice that a flatter EOS corresponds 
to a flatter mass-radius curve. Next we plot figUHl with constant bag pressure of 170MeV, 
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FIG. 16. Gravitational mass (in solar mass) with radius plot of a star sequence with constant 
bag pressure of 170MeV, without magnetic field and with varying magnetic field of field strength 
Bo = lO^'^G and a = 0.005. 

with and without magnetic field. The mass vs radius curve in fig [16] is flatter than fig 
[TSl because this corresponds to the EOS for constant bag pressure, which is much flatter 
than the EOS with varying bag pressure. The varying magnetic field has Bq = lO^^G and 
a = 0.005. Initially, the mass for the star with magnetic field is higher, but the maximum 
mass is lower than the non magnetic case, because the non magnetic EOS is steeper than the 
magnetic counterpart. The stiffness (or flatness) of the pressure vs energy density curve for 
a particular EOS is reflected in the stiffness (or flatness) of the corresponding mass-radius 
curve. 

Next in fig[l7l we plot curves with varying bag constant (ITOMeV) for two different value 
of alpha (0.005 and 0.01), with field strength of i?o = 2 x 10^ ''G. Both the magnetic field and 
bag pressure are density dependent. The magnetic field makes the mass-radius curve flatter. 
As the magnetic field variation becomes higher, increasing the magnetic field strength as we 
go inwards, and thereby making the EOS flat. As the EOS becomes flat the mass-radius 
curve also becomes flat, and the maximum mass decreases. T compare the mass dependence 
on varying magnetic field and varying bag pressure we have plotted curves for two different 
set of curves with varying bag pressure 170 and ISOMeV (fig fTSi) . Each set comprising of two 
curves one without magnetic field and one with magnetic field, of strength Bq = 2 x 10^''G. 
The qualitative nature of the curves remains same due the reasons discussed earlier. As it 
has been pointed out, with varying bag constant and varying magnetic field we can have 
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FIG. 17. Gravitational mass (in solar mass) with radius plot of a star sequence having density 
dependent bag pressure 170MeV. Curves are plotted without magnetic field and with same magnetic 
field, of field strength Bq = 2 x 10^^ G but different a. 

2 
1.8 
1.6 



S 1-2 
1 

0.8 
0.6 

8 9 10 11 12 13 14 15 
radius (km) 

FIG. 18. Gravitational mass (in solar mass) against radius plot of a star sequence with two different 
varying bag pressure of 170MeV and ISOMeV. The curves are plotted without magnetic field and 
with magnetic field, of strength Bq = 2 x 10^''G having same a = 0.01. 

mixed phase EOS with bag pressure of 160MeV. In fig[T9|, we have plotted the mass-radius 
curve for Bg = 160MeV, with {Bq = 2 x lO^^G) and without magnetic field. The magnetic 
field is varying having a of 0.01. The maximum mass for this case is obtained without the 
magnetic field effect and the introduction of the magnetic field makes the curve flatter and 
also reduces the maximum mass. The maximum mass of a mixed hybrid star obtained with 
such mixed phase region is 1.84M0. 

Recently, after the discovery of high-mass pulsar PSR J1614-2230 ^ with mass of about 
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FIG. 19. Gravitational mass (in solar mass) against radius plot of a star sequence with density 
dependent bag pressure of IGOMeV. Two curves are plotted, one without magnetic field and one 
with magnetic field, having strength Bq = 2 x 10^ ''G and a = 0.01. 



1.97Mq, the EOSs describing the interior of a compact star have been put to severe con- 



straint. The [l| typical values of the central density of J1614-2230, for the allowed EOSs is 
in the range 2no - 5no. On the other hand, consideration of the EOS independent analysis 



of 56| sets the upper limit of central density at IOuq. For a constant bag pressure, the 
mass of the HS is about 1.5 solar mass (fig [T6|) . With a varying bag pressure, the maximum 
mass limit can be increased. The maximum mass limit of mixed phase EOS star with the 
above given set of EOS, with strange mass of 160MeV is calculated to be 1.84 solar mass. 
The maximum mass for the mixed hybrid star can be increased to 2.01 57|] solar mass with 
s-quark mass of SOOMeV and varying bag pressure of 150 MeV. Therefore the mass limit 
set by the observation of pulsar PRS J1614-2230 can be maintained by the mixed hybrid 
star having density dependent bag constant. But for this particular choice, the mixed phase 
region is very small. It should be mentioned here that this mass limit is only for this set of 
nuclear and quark matter EOSs. Using very stiff EOS sets (hadronic NL3 and quark quark 
NJL model) the maximum mass limit for the mixed hybrid star can be raised much higher 
as pointed by Lenzi & Lugones 



The main aim of this paper was to show the effect of magnetic field on the mixed phase 
EOS and its effect on the maximum mass of a star. We were also interested in showing 
whether simple EOS (hyperonic nuclear and MIT bag quark) can reach the limit set by PSR 
J1614-2230. The other most interesting fact of this calculation is that the mixed hybrid star 
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has radius corresponding to the maximum mass, quite different from the nuclear and strange 
star. They are not as compact as strange stars and their radius hes between the nuclear and 
strange star. It is also clear from our calculation that, if the magnetic field influence the 
EOS only through the Landau quantization, it has a negative effect on the matter pressure 
thereby making the EOS softer, and the star becomes less massive. 



SUMMARY AND DISCUSSION 

To summarize, we have studied the effect of magnetic field on the nuclear and quark 
matter EOS. We have taken into account Landau quantization effect on the charged particles 
of both the EOS. We have considered relativistic mean field EOS model for the nuclear 
matter EOS. For the quark matter EOS, we have considered simple MIT bag model with 
density dependent bag constant. The nuclear matter EOS is much stiffer than the quark 
matter EOS, and so the effect of magnetic field is much more pronounced in the quark matter. 
The magnetic field due to Landau quantization softens the EOS for both the matter phases 
since the magnetic pressure contributes negatively to the matter pressure. Here we should 
mention that the effect of magnetization of matter is important for strong magnetic fields. 



however it is believed that in NS such magnetization is mall 29| . Therefore in our calculation 
we have neglected such effect. 

Glendenning construction 47(], determines the range of baryon density where both phases 
coexist. At densities below the mixed phase, the system is in the charge neutral hadronic 
phase, and for densities above the mixed phase the system is in charge neutral quark phase. 
We have considered density dependent bag pressure, which has been parametrized according 
to the Gaussian form. We have fixed the lowest value of the bag pressure to be ISOMeV, 
known from the experiments 50|. Accordingly, we have also considered varying magnetic 
field. Observationally, the inferred surface magnetic field of a NS may be as high as lO^^G 
and is believed to increase at the center. As the density decreases with increasing radial 
distance, we have taken the parametrization of the magnetic field as a function of density, 
having maximum field strength at the core. Considering density dependent bag pressure 
and magnetic field, we construct mixed phase EOS following Glendenning construction. 

We find that the effect of magnetic field is insignificant unless the surface field is of the 
order of lO^^G. Such constant magnetic field value has no effect on the nuclear matter EOS 
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and has very little effect on the mixed and quark matter EOS. For a varying magnetic 
field whose surface value is lO^^G but whose central value is of the order of lO^^G, we find 
significant effect on the stiffness of the EOS and also on the the extend of the mixed phase 
region in the EOS. As the bag pressure increases the EOS for the quark phase becomes 
soft, and hence more the effect of magnetic field. At the central region, the bag pressure 
decreases but the magnetic field increases, and so their respective effect on the EOS act in 
the opposite direction. 

The magnetic field increases as we go to much higher densities, and so the boundary 
between the mixed phase and the quark phase changes with increasing field strength. As 
the magnetic field increases, the EOS becomes less stiffer and the phase boundary between 
the mixed and quark phase shifts upwards to the higher density value. Towards the low 
density regime of the curve the effect of magnetic field is less pronounced, as the magnetic 
field strength is less and also the nuclear matter EOS is much stiffer. Therefore the phase 
boundary between the nuclear and mixed phase is less affected. 

The maximum mass limit of mixed phase EOS star is also shown in this paper. We obtain 
a significant mixed phase region with central bag constant of 160MeV having s-quark mass 
of 150MeV. For higher s-quark mass (300MeV) we get a small mixed phase region with bag 
pressure ISOMeV. For such a case we find the maximum mass for a mixed hybrid star with 
the given set of EOS is 2.O1M0. The maximum mass is obtained without magnetic field 
effect and the introduction of the magnetic field reduces the maximum mass. Therefore the 
mass hmit set by the observation of pulsar PRS J1614-2230 is maintained by the mixed star 
with varying bag constant. Our calculation also shows that the mixed hybrid star has radius 
(for the maximum mass) quite different from the neutron or strange star, their radius lying 
between the neutron and strange star. 

Observationally, the surface magnetic field of most of the pulsars are in the ranges of 10^ 
to lO^^G. Such fields have almost no effect in the EOS of matter in those stars. However, 
for magnctars the magnetic fields are very high (~ 10^^ — lO^^G). Flux conservation from 
progenitor stars can give rise to magnetic field of field strength 2 — 3 orders higher. The 
mass-radius relationship for a mixed hybrid star is quite different from the pure neutron 
or strange star, and so it is likely to have different observational characteristics. It is also 
clear that magnetars are different from normal pulsars, as they have lesser mass due to 
flatter EOS. It is to be mentioned here that we have only considered effect from Landau 
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quantization and found that they have significant effect on the mixed phase region once 
it is greater than lO^^G. Here we have not considered the effect from anomalous magnetic 
moment. Anomalous magnetic moment stiffens the EOS, but their effect is significant if 
the magnetic field strength is of the order of lO^^G. For such high magnetic fields the NS 
becomes unstable. Therefore our consideration of only the effect from Landau quantization 
seems alright. As the interiors of the compact stars are hidden from direct observation, we 
have to rely only on the observations coming from their surface. Recent developments has 
been made on measuring accurately the mass of compact stars but a exact measurement 
of their radius is still not possible l| . The knowledge of the radius of a compact stars can 
really give us the hint of the matter components at the star interiors, as we have seen here 
that different EOS provide different mass-radius relationship. 
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